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We investigate the optimal control problem for non-Markovian open, dissipative quantum system. 
Optimal control using Pontryagin maximum principle is specifically derived. The influences of Ohmic 
reservoir with Lorentz-Drude regularization are numerically studied in a two-level system under the 
following three conditions: loo <C lo c , loo ~ lo c or loq lo c , where wo is the characteristic frequency 
of the quantum system of interest, and lo c the cut-off frequency of Ohmic reservoir. The optimal 
control process shows its remarkable influences on the decoherence dynamics. The temperature 
is a key factor in the decoherence dynamics. We analyze the optimal decoherence control in high 
temperature, intermediate temperature, and low temperature reservoirs respectively. It implies that 
designing some engineered reservoirs with the controlled coupling and state of the environment can 
slow down the decoherence rate and delay the decoherence time. Moreover, we compare the non- 
Markovian optimal decoherence control with the Markovian one and find that with non-Markovian 
the engineered artificial reservoirs are better than with the Markovian approximation in controlling 
the open, dissipative quantum system's decoherence. 
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I. INTRODUCTION 

The theory of open quantum systems deals with the 
systems that interact with their surrounding environ- 
ments d i i, i, i, @, 0. Such systems are of great 
interest, and these open quantum systems have been ex- 
tensively studied since the origin of quantum theory Q. 
Despite of the noticeable progresses in the theory, many 
fundamental difficulties still remain. One of the prob- 
lem is decoherence (or loss of coherence) due to the in- 
teractions between system and environment. Recently, 
it received intense considerations in quantum informa- 
tion and quantum computation, where decoherence is re- 
garded as a bottleneck to the construction of quantum 
information processor 0, [HI • The persistence of quan- 
tum coherence is relied on in quantum computer, quan- 
tum cryptography and quantum teleportation. And it is 
also fundamental in understanding the quantum world 
for the interpretation that the emergence of the classical 
world from the quantum world can be seen as a decoher- 
ence process due to the interaction between system and 
environment. 

Various methods have been proposed to reduce this 
unexpected effect, such as the quantum error-correction 
code [!, error-avoiding code [1, ||, minimal deco- 
herence model @, Bang-Bang techniques Q, quantum 
Zeno effect (QZE) [l2[ and decoherence-free subspaces 
(DFS) Unfortunately, all of these schemes cannot 

suppress the unexpect effect successfully for accessorial 
conditions are needed. Altafmi [l4[ pointed out that the 
irreversible decohering dynamics is uncontrollable under 
coherent control. Optimal control technique, which has 
been successfully studied in chemical systems [HI, [H, G3 
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and classical systems [18[, has been exploited to control 
the quantum decoherence @, [H, , where an optimal 
control law was designed to effectively suppress decoher- 
ence effects in Markovian open quantum systems, dy- 
namic coupling in the spin-boson model, and time opti- 
mal control respectively. In this paper, we consider the 
optimal decoherence control problem in non-Markovian 
quantum open system. 

Markovian approximation is used under the assump- 
tion that the correlation time between the systems and 
environments is infinitely short [H, 0, [2l|. For neglect- 
ing the memory effect, the Lindblad master equation has 
been built. However, in some cases, such as quantum 
Brownian motion(QBM) [22| and a two-level atom in- 
teracting with a thermal reservoir with Lorentzian spec- 
tral density [23[, an exactly analytic description of the 
open quantum system dynamic is needed. Especially in 
high-speed communication the characteristic time scales 
become comparable with the reservoir correlation time, 
and in solid state devices memory effects arc typically 
non negligible. So it is necessary to extensively study 
the non-Markovian master equation. We briefly com- 
pare the non-Markovian dynamics (non-Markovian mas- 
ter equation) with Markovian process (Markovian mas- 
ter equation) in Appendix A. For details one can refer to 
Gardiner's book [2l[ or/and Breuer's [l[. 

In this paper the focus will be on the optimal deco- 
herence control of non-Markovian quantum system, par- 
ticularly the simplest system possible, a two-level sys- 
tem governed by the time-convolutionless (TCL) equa- 
tion. We determine control fields which minimize the 
cost functional suppressing the decoherence process by 
applying the Pontryagin maximum principle (PMP) in 
Ohmic reservoir with Lorentz-Drude regularization in the 
following three conditions: luq <C u> c , luq ~ w c , 3> w c , 
where loq is the characteristic frequency of the quantum 
system of interest and lo c the cut-off frequency of Ohmic 
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reservoir. Thus, lu c >C coq implies that the spectrum of 
the reservoir does not completely overlap with the fre- 
quency of the system oscillator and too lu c implies the 
converse case. With this it is possible to engineer differ- 
ent types of artificial reservoir, and couple them to the 
system in a controlled way. We also compare our results 
with no-control system evolution and the optimal con- 
trol of the open system with Markovian approximation. 
The main result of the paper is that decoherence phe- 
nomenon can be successfully suppressed in the u>q <C uj c 
case. Then this explores the coupling of the system to 
engineered reservoirs [U HH , in which the coupling and 
state of the environment are controllable. This may pave 
a newly way to the realization of the first basic elements 
of quantum computers. 

The paper is organized as follows. We first introduce 
quantum decoherence and the quantum master equation 
for driven open quantum systems. In Sec. Ill we formu- 
late the optimal control formalism and deduced PMP 
with a minimum cost functional. Moreover, we con- 
sider the non-Markovian two-level optimal control prob- 
lem. In Sec. IV, we numerically analyze the optimal 
control of decoherence to the two-level system, and an- 
alyze the difference between Markovian optimal control 
and non-Markovian optimal control from both the system 
time evolution and the power spectrum. Conclusions and 
prospective views are given in Sec.V. 



Hi is the control Hamiltonian adjusted by the control pa- 
rameters Ui(t), and Ui(t) represents the control field. The 
Hamiltonian of the environment is assumed to be com- 
posed of harmonic oscillators with natural frequencies u>i 
and masses to,-, 
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where {x\ , x% , •• ■ , xn,Pi,P2, • • • , Pn) are the coordinates 
and their conjugate momenta, and the Planck constant h 
is assigned to be I. The interaction Hamiltonian between 
the system S and the environment B is assumed to be 
bilinear 



Hint = «y]A n ®B n . 



(5) 



The interaction Hamiltonian in the interaction picture 
therefore takes the form 



Hint® = e^+ H ^H mt e-^+H B) t 

= a'£n A n(t)®B n (t), 



(0) 



where 



A n {t) 
B n (t) 



e %Hstj^ g—fflst 
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II. MODELING THE QUANTUM 
DECOHERENCE CONTROL SYSTEM 

Consider a quantum system S embedded in a dis- 
sipative environment B and interacting with a time- 
dependent classical external field, i.e., the control field. 
The total Hamiltonian has the general form 



H to t — H o + Hb + H int 
= H s + H c (t) + H B 



Hi, 



(1) 



where H$ is the Hamiltonian of the system, Hc{t) the 
Hamiltonian of the control field, and Hb the bath and 
Hint their interaction that is responsible for decoherence. 
The operators Hg and Hb act on Hs and Hb, respec- 
tively. The operator Hc(t) contains a time-dependent 
external field to adjust the quantum evolution of the sys- 
tem. One of the central goals of the theoretical treatment 
is then the analysis of the dynamical behavior of the pop- 
ulations and coherences, which are given by the elements 
of the reduced density matrix, defined as 



Ps(t) = tr B [ptot(t)}, 



(2) 



where p to t is the total density matrix for both the system 
and the environment, trs the partial trace taken over 
the environment. The driven model consists of a TV-level 
system interacting with a thermal bath in the presence 
of external control field 0, , and the Hamiltonian is 



H c {t) =^Ui{t)Hi 



(3) 



The effect of the environment on the dynamics of the 
system can be seen as a interplay between the dissipa- 
tion and fluctuation phenomena. And it is the general 
environment that makes the quantum system loss of co- 
herence (decoherence). In general the decoherence can 
be demonstrated as the interaction between the system 
and environment. Then the reduced density matrix of 
the system can evolve into the form 



Pt 



\a n ){a r , 



(7) 



which describes a statistical mixture of nonintcrfcring 
states. Thus, a commonly proposed way to analyze de- 
coherence is by examining how the nondiagonal elements 
of the reduced density matrix evolve under the master 
equation. 

In the present work, we shall concentrate on opti- 
mal control of the decoherence effect in open quantum 
system. The kinetic equation of strong coupling non- 
Markovian quantum system is the following exact time- 
convolutionlcss (TCL) form of the master equation, 

j^pit) =Y d u i £i{t)Vp{t)+1C{t)Vp{t)+l{t)Qp{t ), 

i 

(8) 

with the time-local generator, called the TCL generator 



— JH s t 



H t e~ lHst , K(t) = aVC(t)[l - ^(t)]- 1 ^, 

(9) 
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and the inhomogcncity 

J(t) = aV£(t)[l - mr'git, to)Q, (10) 
where H(t) is the superoperator 

E(t) = a/ dsQ{t,s)QL{s)VG{t,s). 

For detail one see Appendix A and/or [l[ 

In order to facilitate the calculations, we will convert 
the differential equation (8) from the complex density 
matrix representation into the so-called coherent vec- 
tor representation 

[IB EH- Firstly, we choose an or- 
thonormal basis of N x N matrices {(/, ^j)}j=i,2,--- N 2 -i 
with respect to the inner product (X,Y) = tr(X^Y), 
where / is the N-dimensional identity matrix and Clj 
are JVxJV Hermitian traceless matrices. In particular, 
the Hermitian density matrix p can be represented as 
p = jjl + J2i x >, ' ^jj where x = {x\,x 2 , ■ ■ ■x N 2„ 1 ) T is 
a real (N 2 — 1) dimensional vector, called the coherent 
vector of p. This is the well-known Bloch vector repre- 
sentation of quantum systems. Thus the master equation 
(8) can be rewritten as a differential equation of the co- 
herent vector: 

x(t) = O x(t) + ££=i «i(*)0ia;(t) +Lx(t)x(t) + L 2 (t), 

(11) 

with the initial condition, 

x(to) = x 

where Oo, Oi € so(N 2 — I) are the adjoint representation 
matrices of —iHq, —iHi respectively, and xq is the coher- 
ence vector of poi and the term L±(t)x(t) represents the 
decoherence process, k is the number of control fields, 
J"^ UiHi adjusts the quantum evolution such that the 
coherence is conserved. 



III. QUANTUM OPTIMAL CONTROL 
PROBLEM 

A. General Formalism 

As well known, the evolution of the state variable 
x{t) governed by the master equation (fTTj) depends not 

I 



only on the initial state xq but also on the choice of 
the time-dependent control variable u(t). Some earlier 
works to these control problem are listed in the reference 
[30l . I3II |32| . Especially, the exact result was considered 
of the quantum two-state dynamics driven by stationary 
non-Markovian discrete noise in (33| . In this section, we 
are going to suppress the unexpected effect of decoher- 
ence by optimal control technique that wants to force the 
system evolving along some prescribed cohering trajecto- 
ries. The target state chosen is the free evolution of the 
closed system: 



p T (t) = -i[H , p T (t)], (12) 
which is equivalent to x°(t) = e°°(*~*°) xq- The cost func- 
tional is 



J[u(t)} = *[x(t f ),x°(tf)] + I e(x(t),x°(t f ),u(t))dt, 

(13) 



where the functional ^[x(tf), x°(tf)] represents distance 
between the system and objects at final time and the 
functional f* f Q(x(t),x°(tf),u(t)) accounts for the tran- 
sient response with <d(x(t), x°(tf), u(t)) > 0. 

The optimal control problem considered in this paper 
is to minimize the cost functional J[u(t)] with some dy- 
namical constraints. That is, our problem is 



min ueU[t0 , tf] J[u(t)} = *[x(t f ),x°(t f )}+$e(t)dt, 

x(t) = O x(t)+J2i=iUi(t)Oix(t) 

+ Lt(t)x(t) + L 2 (t), 
x(t ) = x , t e [to,tf], 

(14) 

where ( %f[t ,t f ] = {u(-) : [to,tf] — > andu(-) piecewise 
continuous. 

Using the Pontryagin's maximum principle [f8j . the 
optimal solution to this problem is characterized by the 
following Hamilton- Jacobi-Bellman(HJB) Equation 



% + min u ^ [tQtf] {O x(t) + E- = i Ui(t)Oix(t) + Lt^t) + L 2 (t) + <3(x(t),x°(t f ) 7 u(t))} = 0, 
J(x(t f ),t f ) = ^[x(t f )}. 

I 

In general, it is usually difficult to obtain the analytic solution. To illustrate this method and give more insight, 
solution. Nevertheless, one can always have numerical we will consider this problem for the non-Markovian two- 



4 



level system in the following. 



B. Optimal Control of non-Markovian Two-Level 
System 

In this subsection we consider the decoherence of two- 
level system whose controlled Hamiltonian is 

H a = ^{uj a cr z + u x (t)a x + u y (t)a y }, (16) 

where <7fc with k = x,y,z are the Pauli matrices; luq is 
the transition frequency of the two-level system, and u(t) 
is the modulation by the time-dependent external con- 
trol field. In fact, the free Hamiltonian is H$ = \<jJ<j z . 
Then the control Hamiltonian can be described by <x x , cr y 
according to Cartan decomposition of the Lie algebra 
su(2), which was discussed by Zhang et. al in details 
[28l ]. This is the standard model for atom-field interac- 
tion [HiillElll. 

In our two-level system the assumed bilinear interac- 
tion between the system S and the environment B can 
be written as 



Hi n t = a ( a a 



B 



tr_ ® B^) with B = fc t 



(17) 

where a± = (a x ± ia y )/2, the raising and lowering op- 
erator respectively, fcj is the coupling constant between 
the spin coordinate and the zth environment oscillator, 
and a, is the annihilation operator of the ith harmonic 
oscillators of the environment. The coupling constants 
enter the spectral density function J(u>) of the environ- 
ment defined by 



J(cj) 



rriiUi 



-5{lo - uJi) 



(18) 



and the index i labels the different field models of the 
reservoir with frequencies w*. In the continuum limit the 
spectral density has the form 



J( w )=^(-r- 1 exp( ), 



(19) 



where lo c is a cutoff frequency, and rj a dimensionless cou- 
pling constant. The environment is classified as Ohmic, 
sub-Ohmic, and sup-Ohmic accor ding to n = 1, < n < 
1, and n > 1, respectively [iol. l4ll. |43| . 

In this case, the open quantum system can be written 
as follows 0, SI, S3 

PS = -fao[&z,Ps] - %U x (t)[(T x , p S ] ~ %Uy(t)[<T y , ps] 
. A(t)+7(t) ro -i 

_l — ^ \ n > {2a + p S (T- - (J-cr+Ps - Psv-v+l- 

, (20) 

For convenience we map the density matrix of 
the two-level system onto the Bloch vector x(t) = 



(xx(t),x 2 (t),x 3 (t)) T G R 3 defined by x(t) = Tr[ap(t)}, 
which implies that 



xi(t) = poi(*) + Pw(t), 
x 2 (t) = i(poi(t) - Pio(t)), 
X3(t) = poo(t) - Pn(t)- 



(21) 



Then the explicit equations of motion for the compo- 
nents of the Bloch vector read 

x\(t) = -A(t)xi(t) - uj a x 2 (t) + x 3 (t)u y (t), 
x 2 {t) = u) a x\{t) - A(t)x 2 (t) - x 3 (t)u x (t), 
x 3 (t) = -2A(t)x 3 (t) - 2 7 (t) + x 2 (t)u x (t) - Xl (t)uy{t), 

(22) 

where the expressions for the relevant time dependent co- 
efficients, up to the second order in the system-reservoir 
coupling constant, are given by [l], EH 



A(i) = J Q drfc(r) cos(wo T ) 
7(*) = Jo dTp,(r) sin(w T) 



(23) 



with 



k(r) — 2 / °° dujJ(uj) coth.[fku/2kBT} cos(wr), , . 
/i(r) = 2 J °° du>J(uo) sin(wr), 



being the noise and the dissipation kernels, respectively. 
The equation (I22p can be written compactly as 



where 



x(t) =A(t)x(t) +B(t), (25) 



-A(t) -w Uy(t) \ 
4(t) = ( ljo -A(t) -u x (t) 

-u y (t) u x {t) -2A(i) / 



and 





B(t) = ( 

-2 7 (<) 

Let the Ohmic spectral density with a Lorentz-Drudc 
cutoff function, 



T/ \ 2 7o u 2 c 
J( w ) = a;— - 

' J' 1 -J- f , J* 



7T cj; + ur 



(26) 



where 70 is the frequency-independent damping constant 
and usually assumed to be 1. u is the frequency of the 
bath, and lo c is the high-frequency cutoff. For this type of 
spectral density the bath correlations can be determined 
analytically as 



+00 



fc(r) = 4fc B Tco c 2 ^ 



n— — 00 



where v n = 2nnkBT and 

ju(r) = 2fiw 2 e _Uc ' T ' sign r. 



(27) 



(28) 
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Then the analytic expression for the dissipation coeffi- and the closed analytic expression for A(t) is |4S 
cient 7 (i) appearing in the equation (|23p is 

2 2 

7 (t) = ^L[l - e ~™ ot cos(woi) - re"™ ot sin(wot)], 
1 + r z 

(29) 

I 



A(i) = a 2 w Y^{coth(7rr ) - cot(7rr c )e Wct [rcos(w - sin(w t)] + ^ cos(u t)[F(-r c ,t) 

+F(r c ,t) - F(»r ,t) - F(-ir ,t)] - £ siii(t^t)[ a ^g^ ) [(r - *)<3(-r 0> t) (30) 
+(r + i)G(r„, t)] + ^[^(-r CJ t) ~ F(r C: t)})}, 



where ro = LOa/^ksT , r c = w c /27rfcsT, r = uj c /ujq, and temperature limit, we have 
F(:M)= 2 F 1 ( a; ,l,l + ;E ,e- l/lt ), ( 3 h 



A(t) = 2a kT ? {1 " 



— ro;oi 



1 + r a 



G(x,i) = 2 F 1 (2,l + x,2 + x, e - Vlt ) 



(32) 



eosfwoi) sin(o;ot)]}. 

r 

(33) 



In the following we consider the optimal control for- 
2 F 1 (a, b, c, z) is the hypergeometric function and takes ma i ism f our two-level system. For simplicity we define 



the form 

_ . , , ab a(a + 1)6(6+1) , 

2 F 1 (a,b,c,z) = 1 + —Z+ , \ \ s „ ' z 2 + 



the cost functional as: 



E 



n=Q 



lie 2!c(c+l)z 2 

(a)n(6)„ Z n 

(c)„ n! ' 



•/[«(*)] 



[(x(i) - x°(t)) 2 + (9u J (t)u(t)]dt (34) 



where # > is a weighting factor used to achieve a bal- 
ance between the tracking precision and the control con- 
where (a) n is a Pochhammcr symbol. Under the high straints. The corresponding Hamiltonian function is 

I 



H{x{t), u(t),X(t), t) = [(x(t) - x°(t)) 2 + 9u T {i)u{t)\ + \{t) T [A{t)x(t) + B(t)] 

= [( Xl (t) - xl(t)f + (x 2 (t) x° 2 (t)) 2 + (x 3 (t) - x° 3 (t)f + 9(ul(t) + «l(t))] 
+X 1 (t)[-A(t)x 1 (t) - uj x 2 {t) + x 3 (t)u y (t)} + \ 2 (t)[u} xi(t) - A(t)x 2 (t) 
-x 3 (t)u x (t)} + X 3 (t)[-2A(t)x 3 (t) - 2 1 (t) + x 2 (t)u x {t) - Xl {t)u y {t)], 



where A(t) = {X\{t), X 2 {t), \ 3 {t)) T is the so-called La- 
grange multiplier and x°(t) = (xi(t),x 2 (t),x 2 (t)) is the 
target trajectory defined by p = —\ [H ,p]. It is easy 
to see that x°(t) = (x® cos cut — x 2 sin ut,Xi sin cut + 
x® cosWjXg). The optimal solution can be solved by the 
following differential equation with two-sided boundary 
values, 



together with 



an 



a 2 n I = d 2 H(x*(t),u*{t),\{t),t) < 

du 2 I* du 2 — ' 



which implies that 



Ux(t) = 2^{A 2 x 3 - A 3 x 2 }, 
%(*) = 2^{-^3Xi - Aix 3 }. 



(36) 



(37) 



x*(t) = %±=A(t)x(t)+B(t) 
A(t)- dH 
x*(0) 



-™ = -2[x(t)-x (t)]-A(t) T X(t) 

- Xq 

= 



The minimum principle requires the solution of the 
complicated nonlinear equations. When there is one and 
(35) only one solution {x(t), X(t)} it is the required optimal 
solution [l8||. In general, it is difficult to obtain the an- 
alytic solution, if possible existence, to the above opti- 
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mal control problem. So numerical demonstration to this 
problem will be considered in the next section. 



IV. NUMERICAL DEMONSTRATION AND 
DISCUSSIONS 

In this section, we use the formalism of the preced- 
ing section to determine the optimal control of the de- 
coherence. Though the spin-bath models of real systems 
are expected to be more complicated than the two-level 
Hamiltonians considered here, we study the system in 
various aspects to understand the effect of this simple 
system on the decoherence control. 

In our simulations, the system parameters are cho- 
sen as following, x(0) = (^, = ^£=-, — strong coupling 
constant a 2 — 0.01, weighting factor 9 = 1, u>q = 1 as 
the norm unit. Moreover, we regard the temperature as a 
key factor in decoherence process. For high temperature 
fcfiT = 300tJo, intermediate temperature ksT = 3u>o, 
and low temperature fc^T = O.Sluq. Another reservoir 
parameter playing a key role in the dynamics of the sys- 
tem is the ratio r = uj c /^o between the reservoir cutoff 
frequency uo c and the system oscillator frequency ujq. As 
we will see in this section, by varying these two param- 
eters ksT and r = to c j 'u>o, both the time evolution and 
the optimal control of the open system vary prominently 
from Markovian to non-Markovian. 



Figure 1 shows optimal control of decoherence for 
r«l, r = 1, and r > 1 in high temperature reser- 
voir. All of these contain solid line for free evolution, 
dashed line for Markovian optimal control, dotted line 
for non-Markovian optimal control, and dash-dotted line 
for target trajectory. We can see clearly that the decoher- 
ence can be controlled perfectly inrCl reservoir. From 
Figure 2 we can see that the decoherence time Td can 
be delayed for a long time and its amplitude amplified 
heavily with the non-Markovian control. On the other 
hand, Figure 3 shows that the non-Markovian control 
field is changed more rapidly than the Markovian control 
field and the frequency of non-Markovian is more plenty 
than the Markovian, which helps to understand that the 
non-Markovian case is done better than the Markovian 
case and implies that it is necessary to consider the non- 
Markovian case. 

From Figure 1 we can also see that either Markovian or 
non-Markovian optimal control cannot do well when r = 
1 or r — 10. As we discussed before, diffusion is always 
dominant under the high temperature. In the case r«l, 
A(t) > is always true (48|. However, Maniscalco, et. 
al. [48| showed that if r > 0.27 the diffusion coefficient 
A(t) < 0, and the system becomes non-Lindblad. It 
implies that the environment induced fluctuations will 
be large enough. So our control field is negligible when 
comparing with the high-frequency harmonic oscillators 
of the reservoir. 



A. High temperature reservoir 



B. Lower temperature reservoir 



For high reservoir temperature, diffusion coefficient 
A(t) (30) has the approximation form (33), which plays 
a dominant role since A(i) S> 7(i). Note that, for 
time t large enough, the coefficients A(t) and 7(f) can 
be approximated by their Markovian stationary values 
Am = A(t — > oo) and 7m = "f(t — > oo). From eqs.(29) 
and (30) we have 



2 ? 

aw or 



and 



7m = 



Hi = a u — — o coth(7rr ). 
1 + r z 



(38) 



(39) 



Then, under high temperature, noting 

, . , 1 2kT 

COth(7T7"o) ~ 1 H ~ , 

7rr loo 



As decreasing temperature, the amplitude of A(t) be- 
comes smaller and smaller and j(t) becomes larger and 
larger, which is not negligible anymore. There exists a 
time which relate to both the temperature and the ratio 
such that after the time the combination of dissipation 
and diffusion coefficient A(t) — 7(f) < 0, which changes 
the properties of the control system (35). 

Figure 4 shows the non-Markovian optimal control and 
Figure 5 their power spectrum for intermediate temper- 
ature, and Figure 6 and Figure 7 for low temperature. 
At intermediate temperature the non-Markovian optimal 
control plays little role especially in Figure 4(b)and Fig- 
ure 4(c). Note that in Figure 4(a) and 6(a) the free 
evolution is with little decoherence. We note that the 
optimal control does well at low temperature in Figure 
6. In Figure 6, both Markovian and non-Markovian play 
an important role in controlling the decoherence in both 
r = 1 and r = 10. They can make the quantum coherence 
persistence for a long time. 



A H J = 2c?kT-^— 2 . (40) 
1 + r -2 



C. Engineering Reservoirs 



During the last two decades, great advances in laser 
Inserting Eqs.(38) and (40) into Eqs.(35) one can easily cooling and trapping experimental techniques have made 
get the Markovian optimal decoherence control. it possible to trap a single ion and cool it down to very 



7 




FIG. 1: (Color online) Surviving coherence in off-diagonal matrix elements vs time t [Eq.35] under high temperature environment, 
without control action (black solid line), Markovian optimal control(blue dashed line), non-Markovian optimal controlled dotted 
line), and target trajectory(crimson dash-dotted line) at r = 0.1, r = 1, r = 10 respectively. 
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FIG. 2: (Color online) Comparing Markovian optimal control, non-Markovian optimal control with no control under high 
temperature environment, without control action(black solid line), Markovian optimal control(blue dashed line), non-Markovian 
optimal controlled dotted line), and target trajectory(crimson dash-dotted line) at r — 0.1. 



low temperature. These cold trapped ions are the fa- 
vorite candidates for a physical implementation of quan- 
tum computers and realization of the quantum cryptog- 
raphy and quantum tclcportation. All of these rely on the 
persistence of quantum coherence. f24l . [25j are the recent 
experimental procedures for engineering artificial reser- 
voirs. They showed how to couple a properly engineered 
reservoirs with the trapped atomic ion's harmonic mo- 
tion. They measured the decoherence of superpositions 
of coherent states and two-Fock-state superpositions in 
the engineering artificial reservoirs. Several type of en- 
gineering artificial reservoirs are simulated, e.g., a high- 



temperature amplitude reservoir, a zero-temperature am- 
plitude reservoir, and a high-temperature phase reservoir. 

From above discussions we find that our optimal deco- 
herence control fields do well in the engineering artificial 
reservoirs. 

TABLE I. Controllability 



T&r 


Low T 


Med T 


High T 


r=0.1 


slow decay 


slow decay 


controllable (non) 


r=l 


controllable 


uncontrollable 


uncontrollable 


r=10 


controllable 


uncontrollable 


uncontrollable 



Table I shows the controllable property of non- 
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FIG. 3: (Color online) Comparing Markovian optimal control with non-Markovian one in r=0.1 high temperture reservoir. 
Markovian optimal control u^blue solid line) ity (blue dash-dotted line) , non-Markovian optimal control u x (red solid line) u v 
(red dash-dotted line). 



r=0.1 , medium temperature 




r— 1 . medium temperature 



edium temperature 





FIG. 4: (Color online) Surviving coherence in off-diagonal matrix elements vs time t [Eq.35] under medium temperature environ- 
ment, without control action(black solid line), Markovian optimal control(blue dashed line), non-Markovian optimal controlled 
dotted line), and target trajectory(crimson dash-dotted line) at r = 0.1, r = 1, r = 10 respectively. 



Markovian open, dissipative quantum system. When 
r < 1 the system free evolution is with little decoher- 
ence at low and intermediate temperature and our op- 
timal control plays an important role in controlling the 
decoherence phenomenon at high temperature. More- 
over, when r 1 and r = 1 our optimal control also 
plays an important role in controlling the decoherence 
phenomenon at low temperature. They indicate that 
these engineered reservoirs could be designed that the 
coupling and state of the environment can be controlled 
to slow down the decoherence rate and delay decoherence 
time. 



V. CONCLUSIONS 

In the present work, we have studied the optimal con- 
trol of the decoherence for the non-Markovian open quan- 
tum system. In the general formalism we proposed the 
optimal control problem and derived the corresponding 
Hamilton- Jacobi-Bcllman equation. Usually this kind of 
problem is difficult to be analytically solved. Then we 
considered this problem in the non-Markovian two-level 
system. Through transforming its master equation into 
the Bloch vector representation we obtained the corre- 
sponding differential equation with two-sided boundary 
values. 

Finally, we numerically studied the non-Markovian 
decoherence control for three different conditions, i.e., 
luq <C uj c , luq lu c , ujq lu c in the Ohmic environment 
whose spectral density is with a Lorentz-Drude cuttoff 
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FIG. 5: (Color online)non-Markovian optimal controls and their power spectrum in medium temperature reservoir for r = 0.1, 
r = 1, and r = 10 respectively. Non-Markovian optimal control u x (red solid line) u y (red dash-dotted line). 



r=0.1 , low temperature system evolution 




r=1, low temperature system evolution 
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r=10, low temperature system evolution 
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FIG. 6: (Color online) Surviving coherence in off-diagonal matrix elements vs time t [Eq.35] under low temperature environment, 
without control action(black solid line), Markovian optimal control(blue dashed line), non-Markovian optimal controlled dotted 
line), and target trajectory(crimson dash-dotted line) at r = 0.1, r = 1, r = 10 respectively. 
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non— Markovian control 
in low temperature, r=0.1 




control power spectrum 
in low temperature, r=0.1 




non-Markovian control 
in low temperature, r=1 




control power spectrum 
in low temperature, r=1 




non— Markovian control 
in low temperature, r=10 




control power spectrum 
in low temperature, r=10 



0.5 




0.5 



FIG. 7: (Color online) non-Markovian optimal controls and their power spectrum in low temperature reservoir for r = 0.1, 
r = 1, and r = 10 respectively. Non-Markovian optimal control u x (red solid line) u y (red dash-dotted line). 



function. Our numerical results indicated that the de- 
coherence dynamics behaves differently for the different 
environmental condition which leads to significant dis- 
tinctness in the time dependent behavior of the dissipa- 
tion function 7(i) and A(t). We regarded temperature 
as a key factor in the decohcrence effect and showed that 
the decoherence can't be controlled effectively in high 
temperature for both the Markovian and non-Markovian. 
Comparing with the Markovian approximation we be- 
lieved that it is necessary to consider the non-Markovian 
quantum system. Most of all, we analyzed the short time, 
moderate time, and long time decohcrence control behav- 
iors for r = 0.1, which implies uj c <wo. In this case the 
decohcrence can be controlled effectively, which may in- 
dicates that the decoherence rate can be slowed down 
and decoherence time can be delayed through designing 
some engineered reservoirs proposed by Myatt et. al. 
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A. Comparing the non-Markovian dynamics with 
the Markovian dynamics 

1. Quantum Markovian Process and Markovian Master 
Equation 

Quantum Markovian process or Markovian approxima- 
tion is widely used in open quantum system, typically 
in interaction of radiation with matter (weak coupling); 
quantum optics and cavity-QED (weak damping); quan- 
tum decoherence; quantum Brownian motion (high tem- 
peratures); quantum information; quantum error correc- 
tion; stochastic unravelling (Monte Carlo simulations); 
laser cooling (Levy statistics of quantum jumps), and so 
on. The essence of quantum Markovian process contains 
three assumptions: 

• (i)The initial factorization ansatz (Feynman- 
Vernon approximation). At time t — the bath 
B is in thermal equilibrium and uncorrelated with 
the system S 

Ptot{0) = ps{0)<E)p B ; (41) 

• (ii)Weak system-bath interaction (Born approxi- 
mation). 

• (iii)Markovian approximation. The relaxation time 
tb of the heat bath is much shorter than the time 
scale tr (tb <C Tr) over which the state of the 
system varies appreciably. 



11 



Then it induced the dynamical map $ t : 

ps(0) -► Ps(t) = $ t p s (0) = tr B [U t {p s (0) <g> p B )U}}. 

(42) 

With some conditions, like completely positive and Her- 
miticity and trace preservation we get a quantum dynam- 
ical semigroup: $ t = exp[£t], which implies the Marko- 
vian master equation: 



dt 



Pa (t) = Cp s (t), 



(43) 



where generator of time evolution is in Lindblad form: 
Cps{t) = -Uh s ,ps] +^2^[a lP sa\ - \{a\ai, p s }\. 



effective equations (Master equations). As far as we 
known, there are two ways to derive the master equa- 
tion. One is called the path-integral method by Halliwcll 
et.al I27l . Hu et. al [39|, Ford et. al47|, and Karrlein 
et.al[40|, the other is t he p rojection operator method by 
H.P. Breuer 0, 0, |H, Hf . 

The projection operator method is also called the 
Nakajima-Zwanzig projection. The basic idea of the tech- 
nique is to define a map V as 



Vp = tr B {p} ® Pb, 



(45) 



where ps is a fixed environment state and the map V 
is a projection super-operator acting on operators, i.e., 
V 2 = V . Its complementary projection is 



(44) 



Q = l-V, 



(46) 



2. Non-Markovian Dynamics and Non-Markovian Master 
Equation 

Non-Markovian dynamics system is not a new research 
problem, but recently it received considerable considera- 
tion [l], 0, HI, HI, . Comparing with the Markovian 
dynamics it has three properties: (i)Semigroup prop- 
erty violated: slow decay of correlations, strong mem- 
ory effects; (ii)Initial correlations: classically correlated 
or entangled initial states; (iii)Strong couplings and low 
temperatures, with which we can studied the short-time 
behavior and exact evolution of quantum decoherence. 
With the help of these three properties we can derive 

I 



where X is the identity map. Thus the Nakajima-Zwanzig 
equation can be derived [461 ]: 

j t V P {t) = j dsK(t,s)Pp(s)+l(t)Qp(0), (47) 

where K(t,s) is the memory kernel. To second order 
in the coupling constant the general form of the master 
equation can be approximated by 

f t Vp(t) = K(t)Vp(t)+l(t)Qp(0). (48) 
In general, the TCL generator is 



JC(t)p s = - l -[H s {t),p s ]+Y}C l (t)PsD\{t)+D l {t) Ps Cl{t)] - \Y J {D\{t)C i {t) +CjA(*),Ps}, (49) 

i i 

i 

where Cj(i) ^ Di(t), which means that it is not in the Lindblad form. 



[1] H.P. Breuer, and F. Petruccione, The Theory of Open 
Quantum Systems, Oxford University Press, Oxford, 
2002. 

[2] J. Zhang, C.W. Li, R.B. Wu, T.J. Tarn, and X.S. Liu, J. 

Phys. A: Math. Gen. 38, 6587(2005). 
[3] S. Lloyd, and J-J. E. Slotine, Phys. Rev. Lett. 80, 

4088(1998) 

[4] P. Zanardi, and M. Rasetti, Phys. Rev. Lett. 79, 
3306(1997) 

[5] I. L. Chuang, and Y. Yamamoto, Phys. Rev. A. 52, 
3489(1995) 

[6] R. Alicki, M. Horodecki, P. Horodecki, and R. Horodecki, 
Phys. Rev. A, 65, 062101(2002). 



[7] P. Zanardi, Phys. Rev. A 63, 012301(2001). 
[8] J. Von. Neumann, Mathematical Foundations of Quan- 
tum Mechanics, Princeton University Press, 1955. 
[9] M.A. Nielsen, and I.L. Chuang, Quantum Computation 
and Quantum Information, Cambridge University Press, 
Cambridge, 2000. 
[10] M.B. Mensky, Quantum Menseurements and Decoher- 
ence: Models and Phenomenology, Kluwer Academic 
Publishers, Dordrecht, 2000. 
[11] PW. Shor, Phys. Rev. A 52, 2493(1995) 
[12] L. Viola, and S. Lloyd, Phys. Rev. A 58, 2733(1998). 
[13] B. Elattari, and SA. Gurvitz, Phys. Rev. Lett. 84, 
2047(2000). 



12 



C. Altafini, J. Math. Phys. 44, 2357 (2003). 

H. Rabitz, R. de ViVie-Riedle, M. Motzkus, and [33 
K. Kompa, Science 288, 824(2000). [34 
R.J. Levis, G.M. Menkir, and H. Rabitz, Science 292, 
709(2001). [35 
H. Rabitz, Science 299, 525(2003). 
V.F. Krotov, Global Methods in Optimal Control The- [36 
ory, Marcel Dekker INC, New York, (1996). 
H. Jirari, and W. Potz, Phys. Rev. A 74, 022306(2006). [37 

D. Sugny, C. Kontz, and H. R. Jauslin, Phys. Rev. A 76, 
023419(2007). [38 
C.W. Gardiner, and P. Zoller, Quantum Noise (2nd 
Edition), Springer- Verlag, Berlin Heidelberg New York, [39 
2000. 

W.H. Zurek, Rev. Mod. Phys. 75, 715(2003). [40 

B. M. Garraway, Phys. Rev. A 55, 2290(1997). 

C. J. Myatt, B. E. King, Q. A. Turchette, C. A. Sack- [41 
ett, D. Kielpinski, W. M. Rano, C. Monroe, and D. J. 
Wineland, Nature(London) 403, 269(2000). 

Q. A. Turchette, C. J. Myatt, B. E. King, C. A. Sack- [42 
ett, D. Kielpinski, W. M. Rano, C. Monroe, and D. J. [43 
Wineland, Phys. Rev. A 62, 053807(2000). 

A. Shabani and D. A. Lidar, Phys. Rev. A 71, [44 
020101(2005). 

J.J. Halliwell and T. Yu, Phys. Rev. D 53, 2012 (1996). [45 
J. Zhang, R.B. Wu, C.W. Li, T.J. Tarn, and J.W. Wu, 
Phy. Rev. A75, 022324(2007). [46 
K. Blum, Density Matrix Theory and Applications, [47 
Plenum Press, New York, (1981). 

M. Thorwart, L. Hartmann, I. Goychuk, and P. Hanggi, [ 
J. Mod. Opt. 47,(2000). 

M. Grifoni, and P. Hanggi, Phys. Rep. 304, (1998). [49 
J. S. Shao, C. Zerbe, and P. Hanggi, Chem. Phys. 



235,(1998). 

I. Goychuk, and P. Hanggi, Chem. Phys. 324, (2006). 
M. Ban, S.Kitajima, and F. Shibata, Phys. Rev. A 76, 
022307(2007) 

M. O. Scully, and M. S. Zubairy, Quantum Optics, Cam- 
bridge University Press, Cambridge, (1997). 
P. Meystre, Atom Optics, Springer- Verlag, New York, 
(2001). 

C. Anastopoulos, and B. L. Hu, Phys. Rev. A 62, 033821 
(2000). 

S. Shresta, C. Anastopoulos, A. Dragulescu, and 

B. L. Hu, Phys. Rev. A 71, 022109 (2005). 

B.L. Hu, J.P. Paz, and Y. Zhang, Phys. Rev. D 45, 2843 

(1992). 

U. Weiss, Quantum Dissipative System, World Scientific 
Publishing, Singapore (1993). 

A. J. Leggett, S. Chakravarty, A. T. Dorsey, 
M. P. A. Fisher, A. Garg, and W. Zwerger, Rev. Mod. 
Phys. 59, 1(1987). 

J. H. An, W. M. Zhang, Phys. Rev. A 76, 042127 (2007). 
S. Maniscalco, S. Olivares, and M. G. A. Paris, Phys. 
Rev. A 75, 062119 (2007). 

H.P. Breuer, B.Kappler, and F. Petruccione, 
Ann. Phys.(N.Y.)291, 36(2001). 

H.P. Breuer, D. Burgarth, and F. Petruccione, Phys. Rev. 
B70, 045323(2004). 

H.P. Breuer, Phy. Rev. A75, 022103(2007). 

G.W. Ford, and R.F. O'Connel, Phys. Rev. D 64, 

105020(2001). 

S. Maniscalco, J.Piilo, F. Petruccione, and A. Messina 
Phys. Rev. A 70, 032113(2004). 

R. Karrlein and H. Grabert Phys. Rev. E 55, 153(1997). 



